Reconstructing the History of Structure Formation using Redshift Distortions 
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Measuring the statistics of galaxy peculiar velocities using redshift-space distortions is an excellent 
way of probing the history of structure formation. Because galaxies are expected to act as test 
particles within the flow of matter, this method avoids uncertainties due to an unknown galaxy 
density bias. We show that the parameter combination measured by redshift-space distortions, 
j^mass pj-Qyjjjgg g^ good test of dark energy models, even without the knowledge of bias or ag^^^^ 
required to extract / from this measurement (here / is the logarithmic derivative of the linear 
growth rate, and ag^'^'^'^ is the root-mean-square mass fluctuation in spheres with radius 8h~^Mpc). 
We argue that redshift-space distortion measurements will help to determine the physics behind 
the cosmic acceleration, testing whether it is related to dark energy or modified gravity, and will 
provide an opportunity to test possible dark energy clumping or coupling between dark energy and 
dark matter. If we can measure galaxy bias in addition, simultaneous measurement of both the 
overdensity and velocity fields can be used to test the validity of equivalence principle, through the 
continuity equation. 

PACS numbers: 



I. INTRODUCTION 

Observations of the accelerating nature of the Universe 
show that there is fundamental physics at work that we 
do not understand [l|, Q- Many possibilities have been 
postulated including a new form of the vacuum which 
is not present in the contemporary high energy physics, 
or a modification of gravity which would revolutionize 
our understanding of space and time. Discovering which 
mechanism is correct is one of the key challenges for 21st 
century science. 

The detection of acceleration was obtained by us- 
ing supernovae as standard candles, and therefore relies 
on measuring the cosmological geometry. The physical 
process causing the acceleration could also affect struc- 
ture formation, which provides a complementary way 
of distinguishing between models. In particular, mod- 
els in which general relativity is unmodified have dif- 
ferent Large-Scale Structure (hereafter LSS) formation 
timescales compared with Modified Gravity (hereafter 
MG) models [3,0]. Because of this, many previous stud- 
ies have considered testing the signature on LSS forma- 
tion predicted by MG models [g, S 0| ■ Direct observa- 
tions of LSS growth as traced by galaxies are of limited 
value because galaxies are not expected to be simple trac- 
ers of the underlying matter density field, although we 
return to this point later. 

Maps of galaxies where distances are measured from 
spectroscopic redshifts show anisotropic deviations from 
the true galaxy distribution. These differences arise 
because galaxy recession velocities include components 
from both the Hubble flow and peculiar velocities from 
the motions of galaxies in comoving space. Although 
these "redshift-space distortions" are a nuisance when 
trying to reconstruct the true distribution of galaxies, 
they provide a mechanism to measure the build-up of 



structure, which drives these peculiar velocities on large- 
scales. 

In linear theory, and in the absence of bias, a distant 
observer should expect a multiplicative enhancement of 
the overdensity field 6 that is proportional to 1 -I- ffP, 
where / is the logarithmic derivative of the linear growth 
rate, and fj, is the cosine of the angle to the line-of-sight 
p. With a local linear bias, the real-space galaxy den- 
sity field is affected, while the peculiar velocity term is 
not, so the multiplicative factor is changed to 1 -I- /?/i^, 
where /3 = f /b. Because of the /i dependence, this in- 
formation can be extracted from galaxy redshift surveys, 
and a number of methods and their application have pre- 
viously been considered. Analyses have been undertaken 
using the 2degree-Field Galaxy Redshift Survey (2dF- 
GRS) [IQl, measuring redshift-space distortions in both 
the correlation function [111 . Il2[ | and power spectrum af- 
ter decomposing into an orthonormal basis of spherical 
harmonics and spherical Bessel functions [ij] . Using the 
Sloan Digital Sky Survey [lj|, an Eigenmode decompo- 
sition has been per formed to separate real and redshift- 
space effects [ia;ll3- I'^ ^ recent paper, these low redshift 
analyses were extended to z ~ 1 |171| using the VIMOS- 
VLT Deep Survey (VVDS) [ll[l^In addition to mea- 
suring /3 at z = 0.8, this work has pushed exphcitly the 
idea of using large-scale peculiar velocities for constrain- 
ing models of cosmic acceleration. These data, and the 
resulting cosmological constraints, are considered further 
in Section IIII Al In particular we argue that it would 
be better to present results in terms of b{z)asiz) and 
f{z)(7s{z), rather than f3 for local linear bias models. We 
show in Section [IIII that bias- independent constraints on 
/(z)(7g(z) are able to discriminate between some mod- 
els of acceleration as well as /(z), which is commonly 
extracted from (3 by applying an independent (and dif- 
ficult) measurement of bias. Galaxy bias measurements 
tend to have the same fractional error as the redshift- 



spce distortion measurements: for example comparing 
redshift-space distortion results from the 2dFGRS|12] 
with bias measurements from measurements of the 3-pt 
function[36]. 

In addition to a direct measurement of f{z)as{z), 
redshift-space distortion measurements can be used to 
test diverse aspects of LSS, as proposed by Song & 
Koyama f^: geometrical perturbations can be recon- 
structed from the evolution history of peculiar veloci- 
ties. With the assumption of an additional measurement 
of galaxy bias, the continuity equation can be tested, 
and anisotropic stress can be constrained. Those diverse 
tests strengthen our power to constrain theoretical mod- 
els, and are considered further in Section [iVl Before we 
do this, we first review the physics that we hope to test 
using peculiar velocities (Section |TT|, and then consider 
the measurements themselves and what they can directly 
tell us about structure formation (Section IIIip . 



II. LINEAR STRUCTURE FORMATION 



where wx is the equation of state, Spx is the perturbed 
pressure and ax is the anisotropic stress. The continuity 
equation, Eq. ([2]), states conservation of local density. 
The Euler equation, Eq. ([3]), represents the conservation 
of local energy momentum, and describes dynamics of 
perturbed fluids sourced by ^. 

The curvature perturbation $ is constrained to the lo- 
cal inhomogeneity via the Poisson equation. 



fc2<i) = iirGNa^Px 5x 



3aH-— 



(4) 



These equations completely determine the dynamical 
evolution of LSS, within a given expansion history H. 

For models based on general relativity with a standard 
dark energy that does not clump on small scales, the scale 
and time dependence of the evolution of perturbations 
are separable in the matter dominated regime. Both 5x 
and 9x are uniquely determined by the expansion history 
H. Thus we are able to trace the evolution of structure 
formation using observations of either Sx or 9x in these 
models. 



A. Basic equations 

In this section, we briefly review the standard deriva- 
tion of the dynamics of the density fluctuations and their 
associated peculiar velocities in a Friedman universe. We 
will contrast this against non-standard models in later 
sections. In the Newtonian gauge, the perturbed metric 
describing local gravitational instability of the energy- 
momentum density fluid is given by 



ds^ = -(1 + 2*)dt2 + a^{l + 2^)dx^ 



(1) 



where scalar perturbations are dominant over vector or 
tensor perturbations. The Newtonian force ^ sources 
the dynamics of the perturbed fluids, while the curvature 
perturbation $ measures the local energy density fluctu- 
ations. In the epoch of non-relativistic particle domina- 
tion, in standard GR models, $ is identical to — ^, and 
there is no anisotropic stress. 

The spatial variation of density fluctuations is ex- 
pressed by the density contrast dx = Spx{t,x)/px{t) 
where X denotes the specific fluid that we are consider- 
ing. In the approximation of negligible irrotational flow, 
the divergence of the peculiar velocity vx , 6'x = V • vx 
can be used to describe the fluid motion. The dynamics 
of scalar perturbations of the fluid X is well described 
by energy momentum fluctuations dx and 9x, and the 
corresponding metric perturbations ^P and <I>. 

The conservation equation V^T^ = gives the set of 
equations that describe the dynamics of fluid X, 



dSx 
dt 

d9x 
dt 



-[l + wx) iH 

a px 



3HwxSx, (2) 



-H{1 - 3wx)9x - 

fc' f Spx 1 
a \ px l + wx 



dwx /dt 

l + wx 



X 



(3) 



B. Dark energy model without clumping 

In standard dark energy models (hereafter sDE), the 
cosmic expansion is accelerated by introducing a homo- 
geneous dark energy component into the Friedman equa- 
tion, which then predicts an expansion rate H, 



H^ = hI 



n3 



a3 



VL 



DE 



q3(1+io_de) 



(5) 



where '5' denotes baryon, 'c' denotes cold dark matter. 

The incoherence between baryons and CDM perturba- 
tions caused by acoustic waves in the early universe is 
removed at redshifts z > 10 ^20], and can therefore be 
ignored in our analysis. We also assume a negligible cos- 
mological neutrino density, consistent with observations 
(e.g. [2l|). Following these approximations, all mat- 
ter inhomogeneities can be assumed to coherently evolve 
with those in the CDM. Thus we can treat all matter as 
a single fluid denoted by 'm'. Eq. ^ and Eq. ([3]) for 
matter fluid 'to' are. 



d5m 
dt 
d9„i 
"dT 



= 



H9m — 



k\ 

a 



(6) 

(7) 



Since we are assuming that the dark energy is homoge- 
neous on the scales of interest, the metric perturbations 
can be simply related to the total matter-energy fluctu- 
ations. 



k"^^ = AnGiqa^PmSr, 



(8) 



We use ACDM model as an example of sDE with the 
cosmological parameters wd^ch"^) — 0.11, Wbi^bh"^) = 
0.021 and h = 0.72. 



The combination of Eq. ^, the anisotropic stress re- 
lation $ — — \I>, and the coupled equations of Eq. ^ and 
Eq. (O provide a unique solution for the formation of 
cosmological structure. See [22| for a review of methods 
to solve these equations. 



models. If we are able to measure both Sm and 0,„ simul- 
taneously, then the non-trivial anisotropic stress can be 
measured, which also distinguishes DGP and sDE mod- 
els. 



C. Modified Gravity: DGP 

The cosmic acceleration may arise from a modification 
of gravity on cosmological scales as in the DGP model 
Q. In this model, we live on the (3 -I- 1)— dimensional 
brane which is embedded in an infinite Minkowski bulk. 
The weakened gravity at cosmological scales induces 
the cosmic acceleration without introducing dark energy. 
The expansion history of DGP model is determined by 
the usual matter-energy density and the crossover scale 
defined as the ratio of 5-dimensional to 4-dimensional 
Planck mass scales r^ = M^ /SM'jj . 



H SttG 



H' = 



N 



(9) 



As an example of DGP models, we take the same Wc 
and Wb as in our ACDM model discussed in section HTBI 
but use a different h = 0.80 in order to provide a nearly 
identical H{a). 

In general the equations of motion of linear perturba- 
tions in DGP is not closed without solving the dynamic 
equation of propagation through the bulk. But in the 
quasi static limit k/aH ^ 1, where the contribution of 
bulk gradient is negligible, the solution can be derived 
from dynamic equations on the brane. In this regime, 
the perturbed potentials in DGP become modified. 
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1 



3/3. 



DGP 



a^PniSm , (10) 



a'^PrnSrn , (H) 



where 



JDGP 



= 1 - 2rrH 1 



H 



(12) 



The effective Newtonian constant in Poisson equation 
fc^$ = 4:TrGeff{a)a'^pmSm becomes Geff{a) = G'iv(l — 
l/3/3riGp)j and a non-trivial anisotropic stress is intro- 
duced 



$ 
* 



l-3/3_ 



DGP 



1 + SpDGP 



(13) 



while there are no changes in the continuity equation and 
the Euler equation [5, 23, 24]. 

The differences between DGP and sDE models lead 
to distinct evolution of dm and 9m even for models with 
identical expansion histories. Coupling measurements of 
either Sm or 9m with observations of the geometrical evo- 
lution has the potential to distinguish DGP from sDE 



D. Darlc energy model with clumping 

If the dark energy can support long-lived fluctuations, 
then the galaxies will trace the total density fluctua- 
tions St rather than those just in the matter Sm- In 
these clumping dark energy models (hereafter cDE), the 
baryons will fall into the potential wells created by both 
the dark matter and energy. Given the different physical 
behaviour of dark matter and dark energy St will not 
be a simple linear function of Sm ■ The Poisson equation 
with dark energy clumping will be given by 



fc^<i> = AttGno^PtSt ■ 



(14) 



where pxSr = PmSm + PdeSde UM- 

When we have multiple components with different 
equations of state, the peculiar velocities of the differ- 
ent components will not, in general, be the same at any 
spatial location. Consequently, we will measure 9m from 
galaxy redshift-space distortions, rather than 9t- The 
continuity equation holds for matter and DE separately, 
although it does not hold if we mix components. So, for 
example, 



at a 



(15) 



in general. If we can measure St and 9m ^ we can dis- 
tinguish between DGP and cDE by using Eq. (fT5|) . For 
example, if the contribution from dark energy fluctua- 
tions are accidentally identical to the effect of Geff of 
DGP, 



Geff (a) = Gn 



PtSt 



(16) 



and the anisotropic stress causes a similar change to ^ in 
Eq. (fT3|) of DGP, then there is no difference in structure 
formation as measured by the evolution of 9m- However, 
if we observe St, then we can break the consistency of 
Eq. ([T5]), leading to a possible method for distinguishing 
between these models. 



E. Interacting dark energy model 

Current observations allow a coupling between dark 
matter and dark energy [2a, [2^ , although the coupling 
between baryon and dark energy is strongly limited by 
current experimental constraints :2!||. For a model with 
interacting dark matter and dark energy (hereafter IDE), 



the set of background continuity equations are, 



dpb 
dt 
dpc 

dt 



SHpb 
iHp, 



= -ip, 



dpv 
dt 



+ 3H{pv+pi,) = 



# 
' dt 
d(l) 



(17) 



where 7 is an arbitrary couphng constant. A positive 
couphng constant represents decay from dark matter to 
dark energy, and a negative couphng constant represents 
decay from dark energy to dark matter. The equation of 
motion of the scalar field is given by 



di2 ^-^ dt d6 



IPc- 



(18) 



The conservation equation leads to a set of dynamic 
equations for baryons and for dark matter. For baryons, 
there is no signature from the coupling in either the con- 
tinuity equation or the Euler equation. For the dark mat- 
ter, the coupling influences the dynamics. 



d5c Oc 
dt a 
dOr d6 



= 



— * - ay^p^bc 
a 



(19) 

(20) 



The continuity equation still holds because the cre- 
ation/destruction rate of dark matter is proportional to 
its current density, but the coupling between dark matter 
and dark energy changes the Euler equation. 

The interaction between dark matter and dark energy 
modifies the dynamics of dark matter, so two test parti- 
cles of baryon and dark matter placed in the same force 
field will respond differently. There is no change in the 
matter continuity equation but, in general, we will mea- 
sure the velocities of baryonic material, while we observe 
galaxies in the dark matter potential wells, so measure 
matter overdensities. The observed continuity equation 
will then be broken in the presence of coupling between 
dark matter and dark energy [30| . 



III. MEASURING PECULIAR VELOCITIES 

If we map cosmological structure by translating ob- 
served galaxy redshifts to distances assuming that they 
are cosmological in origin then peculiar galaxy veloci- 
ties are misinterpreted, leaving an anisotropic galaxy dis- 
tribution. For pairs of galaxies with large separation, 
the peculiar velocities can tell us about the formation of 
large-scale structure^. On small scales, decoherent pe- 
culiar velocities cause Fingers-of-God (FOG), stretching 
compact structures along the line-of-sight. These distor- 
tions depend on the structure of halos and any cosmo- 
logical information is difficult to distinguish from halo 



properties. The redshift-space power spectrum P^(k) of 
a galaxy redshift survey is commonly modelled (9|,|38| as 

P;^(k) = [Pg,(k) + 2p^P,e,(M) + p'Pe.eM F {k^al{z)p') 

(21) 
where p = k\\/k is the cosine of the angle of the k vec- 
tor to the line of sight. Pg, PsgOg, and Pgg^ are the real 
space auto-power spectra of galaxies and 9g, and the cross 
power spectrum of galaxy-^^ fluctuations, respectively. 
ay and F determine the non-linear velocity distribution 
of galaxies in collapsed structures. It is common to as- 
sume an exponential model for the pairwise peculiar ve- 
locities, so F{k'^al{z)p^) = (1 + fc^cr^(2)/i^)li5 although 
Gaussian models have also been considered [3^ . Even in 
the distant-observer limit, the usefulness of this equation 
is limited because Eq. (|2ip is not physically motivated 
[39|. Additionally, a^ is expected to be a function of 
halo mass and redshift. Consequently, in order to ex- 
ploit the precision available from future surveys, it may 
become necessary to model the exact behaviour using 
simulations [40, |4l|, I44I . In this paper, we assume that 
Eq. (|2ip holds, and that any significant deviations can be 
accurately modelled prior to the new data sets becoming 
available. 

A concern is that the distribution of galaxy pair- 
velocities on large scales might not match the distribu- 
tion of velocities in the matter field. For example, the 
velocity power spectrum for peaks in a density field does 
not match that of the mass [43, |4J|' even if they are 
the same locally. In this paper, we follow the standard 
ansatz that this bias is small (e.g. [7]) and assume that 
9g c^ 9. Current simulations indicate that this might 
cause a 10% systematic error on scales 10-200 /iMpc"^ 
[40I . l4l| , although we are confident that this can be accu- 
rately simulated, so will only affect future measurements 
at a lower level. In the following we therefore drop the 
subscripts g and m from 9. 

We ignore any cosmological information in 
F{k'^(j1{z)p^), and treat this component as a "nuisance 
parameter" to be marginalised over, and concentrate 
on the cosmological information in the linear part of 
Eq. ^. 

For linearly evolving density fields, if wc can write 
down a linear mass conservation equation 9 = —a6m, 
then the growth factor of 9, Dg oc at). Defining / = 
dh\D/d\na = D/{DH), ii 6g is perfectly correlated with 
9 everywhere, then the power spectra in Eq. (I2ip have the 
same shape, and we have that 



Pai^) 



Pg{k)[l + 2p^l3 + p^f3^]F 



m{z) 



(22) 



where (3 = f /b. 

Note that writing Eq. (|22|) in terms of (3 suggests that 
the large-scale redshift-space distortions depend on the 
galaxy bias. We expect the motion of galaxies to locally 
match those of the matter field, so two galaxies at the 
same location would have the same peculiar velocities ir- 
respective of their internal properties: they simply act 
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FIG. 1: The top panel shows the percentage difference in / be- 
tween the sDE and DGP model (long dashed line). The sDE and 
DGP models are described in detail in Section IIIICI The back- 
ground expansion has been matched between these two models. 
The bottom panel shows the percentage difference of /ct™'"^'^ be- 
tween sDE and DGP. Here, the long-dashed curve includes CMB 
data, normalising the models at the epoch of last scattering (using 
A? ) , while the dashed curve shows the model normalised using 
a low redshift measurement of cr'^'^^{z = 0) = 0.82, matching the 
5-year WMAP best-fit ACDM value [H. The blue and black er- 
ror bars are estimated from BOSS and EUCLID respectively (see 
Section IlIIB I for details). 



as test particles within the matter flow. This assump- 
tion is part of Eq. ((2T|) : the dependence on bias simply 
comes from expecting a multiplicative correction to the 
observed galaxy power spectrum for the redshift-space 
distortions. If we instead consider modelling an additive 
contribution, then the dependence on bias can be broken. 

Comparing Eq. (f2T|) and (j22|) shows that, although 
many observational studies present results in terms of 
/3, the fundamental constraint from the fi dependence of 
the normalisation of Eq. ((2T|l is on the normalisation of 
Pgo, and the cross power Pgg. If we assume that the con- 
tinuity equation holds, and we have a local linear bias 



so S„ 



h6„ 



then the normalisation of Pgg depends 



on {haf'^^'f, Pgg depends on (/crg"''"")(6CTf''''')> and Pgg 
depends on {fa^'^^)'^- If the redshift-space distortions 
are modeled as a multiplicative component, then we con- 
strain (bag, (3), while if they are modeled as an additive 
component we measure (bag^fas). These constraints will 
be a simple transform of each other, and should give the 
same likelihood surface if we transform into the same ba- 
sis (this will be considered further in future work|54|). 
Previous analyses have, in general, converted observa- 
tional results into a constraint on / using a separate 
measurement of bias (e.g. [12|)- A perfect estimate of 
bias would allow us to break fcr^'^^ and ba^'^^^ to pro- 



FIG. 2: The top panel shows the time evolution of /: solid curve 
is for sDE model, and a long-dash curve is for DGP. The sDE and 
DGP models are described in detail in Section llll CI The bottom 
panel shows the evolution of fcr™'^^^ : the solid curve shows the sDE 
model, the long-dash curve DGP normalized using A? , and the 
short-dashed line, DGP o-™='==(2 = 0) = 0.82. The blue" and black 
error bars are estimated from BOSS and EUCLID respectively (see 
Section llllBl for details). 



vide separate constraints on a™'^^, f and b. Note that /, 
^mass g^j-^j ^ g^j.g gj^ Tcdshift dependent functions. 

In Fig. [T] and Fig. [21 we show that the combined quan- 
tity, fa^'^^^ is sufficient to distinguish between DGP and 
sDE models described in Section IllICi if the fluctua- 
tions are additionally normalised using CMB data. If 
however, we normalise the model predictions for fa'^^^^ 
using a present day measurement of a'^^^^, then the pre- 
dicted differences in / and fa™'^^ between the DGP and 
sDE models are very similar, a consequence of the weak 
difference in the low-redshift evolution of a™'^^ . Normal- 
ising the model prediction for fa^^'^^^ at low redshift is 
a fairer comparison between using the two ways of de- 
scribing the redshift-space distortion constraints, as this 
effectively only uses the ratios between measurements of 
/cr™'"'''' at different redshifts. 

Given that no technique is significantly better at dis- 
tinguishing these models in terms of the statistical error, 
it is worth considering the potential for including sys- 
tematic errors using both techniques. Here, fa^'^^^ wins, 
since it can be measured without knowing the bias b or 
the amplitude of the matter fiuctuations a™'^^'^. 

In this analysis we have compared two models of accel- 
eration. Clearly, future measurements will be used to test 
a wide range of models. It is clear that for some of these 
/ will be superior, while for others using fa'^^^^^ directly 
would provide tighter constraints. In particular, using 
j^mass (iii-ectly can be advantageous for a model such as 
the DGP model tested here in which growth of structure 



formation is suppressed compared with ACDM. This is 
not true for other theoretical models in which growth of 
structure is less suppressed, e.g. normal branch DGP, 
compared with the self-accelerating DGP model consid- 
ered here 1551. 



A. current constraints 

Redshift-space distortions have recently been mea- 
sured at z = 0.77 using the VIMOS-VLT Deep Survey 
(VVDS)[ll. They use crg"^""(2 = 0) = 0.78 ± 0.03 mea- 
sured by the WMAP experiment, extrapolated to the 
survey redshift z = 0.77 cr^^""(z = 0.77) = 0.54 ± 0.03 
(assuming a fiducial ^m — 0.25, f^A = 0.75 cosmology) to 
measure the bias. Using this estimate, they convert their 
measured value of /3 to provide a constraint on /(z). If 
we do not do this bias modelling, then we can reconstruct 
the constraint on fa^^^ from their quoted results. To 
do this, we take their /(z ~ 0.77) = 0.91 ± 0.36 measure- 
ment and multiply by the extrapolated WMAP result 
^mass^^ = 0.77) = 0.54 Doing this, their redshift-space 
distortion measurement becomes fa^'^^ = 0.49 ±0.18 at 
z = 0.77. We have conservatively assumed that all of 
the quoted error comes from the LSS observations rather 
than the extrapolated CMB normalisation and subse- 
quent bias estimate. 

At lower redshift, we use results from a spherical har- 
monics analysis of the 2dFGRS[lj|. This work did not 
provide a direct constraint on fa^^^ at the fiducial red- 
shift of the 2dFGRS, but extrapolated this to z — 0. 
Without loss of generality, we can undo this extrapola- 
tion using their fiducial cosmology. For this cosmology 
{n^ = 0.3, flA = 0.7), D{z = 0.17)/L'(z = 0) = 0.916, 
and /(z = 0.17)//(z = 0) = 1.2. Their extrapolated 
z = measurement of f(J^'^^^ = 0.46 ± 0.06 should there- 
fore be considered a measurement fa'^'^^^ = 0.51 ± 0.06 
at z = 0.17. 

We also include a measurement of redshift-space dis- 
tortions from the SDSS LRG catalogue [16J. This work 
measured /3 = 0.309 ± 0.035 at z = 0.35. To convert 
to a constraint on fa^'^^^, we multiply by ba^'^^^ = 1.43 
from the joint SDSS-f CMB fit presented. Here we as- 
sume that, because of the wide survey geometry in the 
SDSS, the measurement error on (3 is significantly larger 
than that on ba'^'^^^ . This gives a measurement fa^'^^^ — 
0.44 ±0.05. 



B. future constraints 

In order to see how these constraints will improve with 
the next generation of experiments, we have calculated 
expected errors for three types of future experiment. We 
have considered a ground based redshift-survey of galax- 
ies out to z = 0.8, based on the BOSS experiment, part 
of the SDSS-III (see www.sdss3.org). For this survey 
we have assumed that redshifts will be measured for 



1.5 X 10^ galaxies over lOOOOdeg with approximately 
constant number density out to z = 0.7. We also consider 
a next generation space based mission that could mea- 
sure 0.5 X 10^ galaxies in a volume V = Ix 10^^ h~^ Mpc^ 
out to z = 2. Such a survey is proposed has been pro- 
posed to ESA mission as part of the EUCLID mission, 
and is the result of the merging between SPACE and 
DUNE [45, 46,, 47] . We have also considered the type of 
experiment that could be accomplished using a wide-field 
multi-object spectrograph on a 8m-class telescope. Here, 
we took the "fiducial survey" parameters presented for 
the proposed WFMOS instrument Q. These are two 
surveys with number density 5 x 10"'' h'^ Mpc~ , one of 
2 000 000 galaxies with 0.5 < z < 1.3, and one of 600 000 
galaxies with 2.3 < z < 3.3. The numbers and the 
derived expected measurements for all of these surveys 
should provide an approximate guide to the improve- 
ments that are expected beyond current observations. 

To calculate expected errors for these surveys, we use 
the fitting formula presented in [13], which was derived 
from numerical simulations. This formula was presented 
as providing the error on /, and relies on having an inde- 
pendent accurate measurement of a™'^^ (which is itself 
dependent on cosmology through the growth rate, rather 
than its derivative). In the above discussion we argued 
that it is more consistent to measure and test a constraint 
on fa^^^^^^, and we translate the ^ITJ formula as giving a 
fractional error 
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where V is the volume of the survey (or part of the 
survey under consideration) in (ft.~^Mpc)'^, and Ug is 
the galaxy number density in (/iMpc~^)'^. By perform- 
ing this translation we have conservatively assumed that 
the error predicted on / by [l7| was dominated by the 
measurement of /ct™'*'^'', rather than the translation to 
/. The [17] formula has behaviour very close to Poisson 
with ex l/A^gai, consistent with the idea that we are not 
concerned with the scales on which we see galaxy pairs, 
but simply want amplitude information so all pairs count 
equally, and also assumes that we can recover information 
from both Pgo and Pgg [53]. This formulae is too sim- 
plistic to capture many of the dependencies, such as on 
galaxy bias and power spectrum shape, so we have per- 
formed a full Fisher matrix calculation. This validates 
this simple [l'^ formulae for reasonable values of galaxy 
bias, and galaxy power spectra. As our aim in this paper 
is to present the case for redshift-space distortions as a 
probe of dark energy models, the ,1X] formula is adequate 
for our purpose. We convert to calculate errors on / for 
Figs. [1] and O assuming that 
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the ACDM constraint 
, and assume that the 




with erg = 0.82 ± 0.03, following 
from the 5-year WMAP data [^ 



same level of accuracy can be achieved at low redshift by 
the next generation of lensing experiments. 



C. model comparison 

Current and future constraints on fcr™'^^ are com- 
pared with various models in Fig. [31 The models are 
all normalised to match the amplitude of fluctuations 
seen in the CMB by the 5-year WMAP experiment [STf . 
We see that future experiments will provide an order- 
of-magnitude improvement over current studies, and will 
provide constraints over a range of redshift from a single 
survey. 

When comparing against models, we have not 
marginalised over parameters that are determined by the 
background, but instead fix parameters close to our sDE 
fiducial model in fitting to current geometrical constraint 
given by combination of WMAP, SN and Hubble key 
project. Fig. [3] therefore shows the differences between 
structure growth in the different models, rather than per- 
forming a full Likelihood test given current data. It shows 
what we can tell about structure formation if we have 
perfect information about the background expansion. 

In detail, for DGP, we fix the distance to recombina- 
tion at z\ss = 1088^^2 through the measurement of the 
acoustic peak scale Ia = 302^"'^ and its length calibra- 
tion through the matter density ^mh^ = 0.128 ± 0.008, 
and we fit the DGP cosmological parameter to the com- 
bination of SN and Hq from Hubble Key Project. Open 
curvature is introduced to enhance the fit to data set. 
For cDE, we use the effective dark energy which makes 
the cDE geometrical factor exactly identical to the open 
DGP model used here. For IDE, we fix ^mh^ at its best 
value measured by WMAP, and we allow Hq to be varied 
within 1-CT range of measured Hq = 72 ± 8km/sec/Mpc. 
Here we choose Hq at the edge of upper bound, Hq — 
80km/sec/Mpc, at 7 = 0.3. If 7 becomes bigger than 
this value, then it is expected that it deviates significantly 
from the geometrical measurement as well as structural 
departure considered in this section. 

In this paper we have followed standard convention 
and used a"^^^^ to normalise the matter power spectrum. 
The redshift-space distortion measurements actually con- 
strain fluctuations on a range of scales larger than those 
probed by ctsi but this does not matter if we have a 
constant power spectrum shape. If the power spectrum 
shape is not constant, then the statistic used would have 
to be revised, and the full modelling of P{k) included. 
This is beyond the scope of the current paper, and we 
assume a constant P{k) shape. Here, we have assumed 
that the power spectrum of the initial fluctuations in the 
comoving gauge as measured by WMAP, has amplitude 
A? . = 2.4 X 10"^. We assume that the time-dependent 
growth function of <& is normalized to unity at the on- 
set of matter domination and that the scale dependent 
transfer function is normalized to unity in the fc — > 
limit. For DGP and cDE models, we tune the acoustic 



peak structure to be identical to sDE with same LOm and 
ujb and with fixing angular diameter distance to last scat- 
tering surface. For IDE, we do the same with the extra 
assumption that the earlier DE decay is negligible. This 
is a good approximation for small coupling limit used in 
this paper. The full treatment of the CMB in IDE models 
is beyond the scope of the current paper. 



IV. FURTHER APPLICATIONS OF PECULIAR 
VELOCITY MEASUREMENTS 

We now discuss additional cosmological applications 
that can be developed from peculiar velocity measure- 
ments. This is a development of previous work that con- 
sidered constraints from observations of projected power 
spectra 18] . Here we assume that we have measurements 
of the growth of Pgg. We have a similar measurement 
from the combination of Pgg and Pgg, but as this uses 
the continuity equation for baryons we treat this sepa- 
rately. First, we consider how a redshift dependent mea- 
surement of the amplitude of Pgg can be used to recon- 
struct ^. We quantify this normalisation using cr|, de- 
fined as the rms fluctuations of 6 averaged over spheres 
of radius ?>h~^ Mpc. 

If we are able to probe the history of a"^^^^ , in addi- 
tion to tJg , then we can complete the structure formation 
test. Measuring the derivative of a^'^^^ with respect to 
time would enable us to test the continuity equation by 
comparing it with cTg. This test eliminates many dif- 
ferent theoretical models explaining cosmic acceleration, 
and is potentially more powerful than simply measuring 
only one of these quantities. Finally, the reconstructed ^ 
created from erg can be compared with $ estimated from 
^mass^ which wiU constrain the anisotropic stress. 

Measuring the evolution of a"^"^^ is problematic be- 
cause we cannot observe the mass directly. We could try 
to predict how galaxies trace the matter (e.g. [3l|, [32, 
[33l . [34l . [35j). or could measure this bias, for example us- 
ing higher order statistics (e.g. [3a]). Alternatively, we 
could use weak lensing measurements [3, [4^, [50] ■ We ex- 
pect that these methods will give constraints on bias at 
the percent level, by the time that redshift-space distor- 
tion measurements are obtained from the EUCLID ex- 
periment. Note that the weak- lensing constraints could 
themselves come from an imaging component of the EU- 
CLID mission. For our proposed tests that require bias 
measurements, we therefore assume an expected uncer- 
tainty of A6/6 ^ 0.02. The models used match those 
described in Section fill CI 



A. Reconstruction of perturbed potential 

If (j| can be measured with sufficient accuracy that its 
derivative can also be determined, then we can estimate 
^ from Eq. ([3]). The measured derivative of o'g(z) as a 
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FIG. 3: The solid curve represents sDE, long-dash curve represents DGP, and dotted curve represents IDE. These models 
predict approximately the same background expansion, which is normalised at high redshift to match the fluctuations observed 
in the CMB. The left panel shows the current constraints discussed in the text. The three other panels show simulated data 
for BOSS, WFMOS and EUCLID-type experiments (see text for details). 



function of redshift can be written 
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where i denotes each redshift bin. 
construct an estimator of ^* from 



From this, we can 
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We show estimates of how well future experiments will 
be able to reconstruct ^ using this method in the top 
panel of Fig. 21 If there is no interaction between dark 
matter and dark energy, then (^) = \1/, since there is 
no change from the standard Euler equation. This offers 
a new way to probe the geometrical perturbation other 
than weak lensing experiment. Even with weak lensing 
measurements, it is interesting to have ^ separately, be- 
cause weak lensing probes the combination $ — \l/. 



B. Test on continuity equation 

If we can measure the evolution of a^'^'^ and erg simul- 
taneously, then the continuity equation can be tested. 
The measured a'f^^ at each redshift bin leads us to es- 
timate its derivative in terms of time, which is defined 
by, 
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H, 



(27) 



where i denotes the number of redshift bin. In the middle 
panel of Fig. [H we show the estimated errors in constrain- 
ing the continuity equation with EUCLID. The depar- 
ture from the continuity equation of cDE model having 
the identical structure formation to DGP model (dash 
curve) is detectable. It means that MG models can be 
distinguishable from DE type model, i.e. modified grav- 
ity is detectable. The dotted curve representing IDE has 



detectable departure from the continuity equation. This 
test will be a crucial future test for exotic DE models, 
but requires either galaxy bias to be fully understood or 
modelled. An alternative test using weak-lensing obser- 
vations is possible. 



C. Constraint on anisotropic stress 

With the same assumption as in the previous subsec- 
tion of an independent measurement on bias, we can ex- 
tend the analysis to constrain anisotropic stress. The 
curvature perturbations are determined by 



A:2$« = ATrGNO^PmSl 



(28) 



so a measurement of the normalisation of the density field 
leads to a constraint on $'. Using this measurement 
and the measurement in Section IIV Al we are able to 
constrain anisotropic stress by comparing reconstructed 
^' and $'. We show in the bottom panel of Fig. [H that 
we are able to constrain the level of anisotropic stress 
predicted from DGP model. 



V. CONCLUSIONS 

One of the best ways of measuring this structure 
growth is to use observations of redshift-space distor- 
tions. This is an interesting and timely subject given 
current observations and those planned on a 10-20 year 
timescale [13, [5l|, [5z]. In this paper we have reviewed 
the importance of redshift-space distortion measurements 
given that they provide a measurement of structure 
growth that is independent of galaxy density bias. We 
have argued that peculiar velocity measurements are best 
presented in terms of CTg, or fa^'^^ for models where 
the continuity equation holds. The independence from 
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FIG. 4: Constraints on physical measurements that can be de- 
rived from a EUCLID type survey. The solid lines, long-dash lines, 
dashed lines and dotted lines are sDE models, DGP, cDE model and 
IDE model respectively. The top panel shows the reconstruction of 
^ from (7g , the middle panel shows the test of the continuity equa- 
tion [30| , and the bottom panel shows the constraint on anisotropic 
stress. The errors shown in this figure assume that we can recover 
o-g from both Pgg and Pgg. 



galaxy density bias has not been widely covered in pre- 
vious literature. Most previous analyses have consid- 
ered measuring (3 and ba^^^^, although there are some 
exceptions |l3J. Although extremely simple, we have fo- 
cused on the density bias independent constraint result- 
ing from multiplying these together. The physical origin 
of such a constraint is that linear velocities, which scale 
with the derivative of the growth factor, depend only on 
the matter velocity field. 

The primary conclusion from our work is that con- 



straints on CTg or /cr™^****, are extremely good at helping 
to distinguish between the dark energy models that we 
reviewed in SectionlTTl In fact, as we show in Fig.[TJ these 
constraints are equivalent to similar percent measure- 
ments of / for some models of cosmic acceleration. They 
also have the simplicity of not having to model galaxy 
bias. The simple formula of Tr| has been adapted to 
determine constraints on /cr^^^'^^and shows future con- 
straints in Fig.[3l Although the [I4I formula is simplistic, 
and might not be believed at the percent level, it shows 
that we can expect a huge step forwards in redshift-space 
distortion measurements with the next generation of sur- 
veys. 

Going beyond simply obtaining a single measurement 
of cTg or /cr™''*'^, we have considered how the underly- 
ing perturbation evolution can be tested using peculiar 
velocity measurements. Peculiar velocity measurements 
are important because they can be used to reconstruct 
Newtonian potential '^ which sources the dynamics of 
a galaxy given by Euler equation. Weak-lensing only 
measures ^ in the combination <I> — ^, so redshift-space 
distortions offer a complementary test of perturbations. 

We have considered how peculiar velocities can be used 
to test the continuity equation, which is worthwhile since 
there are many theoretical models which fail to satisfy 
this relation. If dark energy couples to matter, then cur- 
rent constraints show that it must couple to the dark 
matter and not to baryonic material. The coupling of 
dark energy to dark matter modifies the Euler equation 
for dark matter, and breaks the equivalence principle be- 
tween dark matter and baryon. This difference in free- 
fall breaks the continuity equation in which the peculiar 
velocity of matter is estimated using baryons, while we 
consider the growth of fluctuations in all matter. In ad- 
dition, this test can tell if there is dark energy clustering 
which deepens the curvature potential well because we 
measure the peculiar velocity of the matter not of the 
total energy density. Finally we have considered how we 
can constrain the anisotropic stress by comparing $ and 
^, reconstructed from the density fields and peculiar ve- 
locity respectively. 
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